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ABSTRACT 

We report the detection of the lens cluster of the gravitational lens (GL) system Q0957+561 
from a deep observation with the Advanced CCD Imaging Spectrometer on-board the Chandra 
X-ray Observatory. Intracluster X-ray emission is found to be centered 4'f3lJ'3 ^'-^-o 6 

north of image B, nearer than previous estimates. Its spectrum can be modeled well with a 
thermal plasma model consistent with the emission originating from a cluster at a redshift of 
0.36. Our best-fit estimates of the cluster temperature of Tg = 2.09^0^54 keV (90% confidence) 
and mass distribution of the cluster are used to derive the convergence parameter k, the ratio 
of the cluster surface mass density to the critical density required for lensing. We estimate 
the convergence parameter at the location of the lensed images A and B to be ka = 0.22^QQy 
and Kb = 0.211;q;q7, respectively (90% confidence levels). The observed cluster center, mass 
distribution and convergence parameter k provide additional constraints to lens models of this 
system. Our new results break a mass-sheet degeneracy in GL models of this system and 
provide better constraints of ~ 29% (90% confidence levels) on the Hubble constant. We also 
present results from the detection of the most distant X-ray jet (z = 1.41) detected to date. 
The jet extends approximately 8" NE of image A and three knots are resolved along the X-ray 
jet with flux densities decreasing with distance from the core. The observed radio and optical 
flux densities of the knots are fitted well with a synchrotron model and the X-ray emission is 
modeled well with inverse Compton scattering of Cosmic Microwave Background photons by 
synchrotron-emitting electrons in the jet. 

Subject headings: gravitational lensing — galaxies: clusters individual (Q0957-I-561) — X-rays: 
galaxies 



1. INTRODUCTION 

In an earlier paper based on deep ROSAT obser- 
vations, we had reported the 3 a detection of X- 
rays from the cluster of galaxies of the gravitational 
lens system Q0957-F561 (Chartas et al. 1998). This 
cluster contributes to the lensing of a distant z = 
1.41 radio loud quasar. The lensed quasar appears 
as two images separated by Q''17 and denoted as A 
(north) and B (south). An accurate determination 
of the mass distribution of the cluster is essential 
in reducing the uncertainty of the Hubble constant 
as derived from the application of Refsdal's lensing 
method (Refsdal, 1964a, 1964b). Due to the limited 
spatial and spectral resolution and low signal-to-noise 
ratio of the previous X-ray observations, only weak 
constraints could be placed on the cluster properties. 
In particular, we obtained estimates for the conver- 



gence parameter, k, the ratio of the projected two- 
dimensional (2D) surface mass density of the cluster 
to the critical surface mass density of the lens cluster. 
We found k to range between 0.07 and 0.21, when as- 
suming that the cluster center was located 24" away 
from image B, a separation based on optical observa- 
tions of the galaxy members of the cluster (Angonin- 
Willaime et al. 1994). Other methods for estimating 
the convergence parameter k rely on weak-lensing of 
background galaxies (Fischer et al. 1997), measure- 
ments of the velocity dispersion of the lens galaxy Gl 
(Falco et al., 1997; Tonry & Franx 1999; Romanowsky 
& Kochanek 1999) and measurements of the velocity 
dispersion of the lens cluster (Garrett, Walsh, & Car- 
swell 1992, Angonin-Willaime, Soucail, & Vanderri- 
est 1994). The present constraints on k are quite 
poor and have probably underestimated the system- 
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atic errors due to the uncertainty of the location of 
the cluster center and of the mass profile. Barkana 
ct al. (1999) found a strong dependence of k on the 
assumed cluster mass profile and on the distance of 
the cluster center from the lens galaxy. 

A measurement of the convergence parameter k at 
the location of the lensed images with respect to the 
cluster center is needed to break the Hubble constant 
mass-sheet degeneracy (Falco et al. 1985). More 
direct methods of determining the location of the 
center of mass of the cluster are provided by mea- 
surements of the intracluster medium emission in the 
X-ray band, the weak lensing method, and optical 
measurements of the member galaxies. X-ray mea- 
surements with ROSAT did not provide any useful 
constraints on the location of the cluster center. The 
weak lensing method, which determines the mass dis- 
tribution of the lens cluster from the gravitational dis- 
tortion of the images of background galaxies, yields 
a center of mass of the cluster to be located ~ 18" 
east and 13" north of Gl, with a la uncertainty of 
about 15". Angonin-Willaime et al. (1994) derived a 
center for the cluster at 13"7 west, 19"6 south of im- 
age B by proportionally weighting galaxies by their R 
band luminosity. This method was found to be very 
sensitive to the number of galaxies included in the 
calculation. For example, Garrett et al. (1992) using 
a subset of the galaxies had found the cluster center 
at 3"5 west and 1"2 north of Gl. Another method for 
estimating the cluster center is based on the "best- 
fit" values for the lens models that describe this sys- 
tem. In a recent analysis Chae et al. (1999) model 
the cluster contribution by a power-law sphere with 
an extended core. They found a distance between 
the cluster center and the Gl galaxy of 9"^| and a 
position angle 9 (North through East) of ~ 52°. In 
an independent analysis Barkana et al. (1999) model 
the cluster as a singular isothermal sphere and al- 
low the center position and velocity dispersion of the 
cluster to be free parameters. They found a value 
for the distance of the cluster center from the Gl 
galaxy of 13"7 east, 6"9 north and a value for the 
velocity dispersion of 439 km s~^. The large system- 
atic differences between independent estimates of the 
center position of the cluster obtained through mod- 
eling imply that the mass distribution of the lens is 
not accounted for correctly in these models. A recent 
discovery of the lensed images of the host galaxy of 
the quasar Q0957-F561 by Keeton et al. (2000) has 
provided tighter constraints on the global structure of 
the lensing potential. Gravitational-lens (GL) mod- 
els of Q0957+561 developed prior to the discovery of 
the host galaxy failed to reproduce the correct shape 
for the host galaxy arcs and therefore any constraints 
they provide on the value of the center location of the 



cluster are unreliable. Lens models by Keeton et al. 
(2000) that incorporate the constraints provided by 
the host galaxy arcs imply that the cluster potential 
is approximately centered on Gl. 

In this paper we present results from the recent 
observation of the GL system Q0957-I-561 with the 
Chandra X-ray Observatory. Section 2 contains a de- 
scription of the observation and the data reduction. 
In section 3 we perform a detailed spatial analysis 
of the X-ray image of the lens cluster and present a 
more direct measurement of the location of the clus- 
ter center and the cluster core radius. In section 3 
we also briefly describe the detection of an X-ray jet 
corresponding to the radio jet from image A. The 
spectral analyses of the lens cluster, lensed images, 
and X-ray jet are presented in section 4. We revisit 
the determination of the convergence parameter, n, in 
section 5. In section 6 we discuss the implications of 
the Chandra constraints on the Hubble constant and 
conclude with a summary of the results obtained from 
the Chandra observation. We use Hq = 75 km s"^ 
Mpc~^, qq = 0.5, and A = 0, unless mentioned other- 
wise. At the redshifts of 0.36(lens) and 1.41(quasar) 
an angular size of 1" corresponds to length scales of 
~ 4.1 kpc and ~ 5.7 kpc, respectively. 

2. OBSERVATIONS AND ANALYSIS 

Q0957-f561 was observed with the ACIS instru- 
ment (Garmire et al. 2001, in preparation) onboard 
the Chandra X-ray Observatory in 2000 April 4 
for 47,660s. The data were collected on the back- 
illuminated S3 CCD of ACIS. The instrument was 
operated in sub-array mode where a restricted region 
with 254 of 1024 rows of the chip collects data to 
reduce the frame-time from the nominal 3.24 s (full- 
frame) to 0.741 s. The telescope pointing was set ap- 
proximately 1 arcmin off axis. We used the sub-array 
mode and off-axis pointing to reduce the probability 
of multiple X-rays landing within a few CCD pix- 
els during the same CCD frame readout. This effect, 
commonly referred to as pile-up, may lead to spectral 
and spatial distortions and loss of detected events. 
The detected count rates for images A and B are 
0.21 counts s^^ and 0.16 counts s~^, with respective 
estimated pile- up fractions of about 11% and 7% for 
the 0.741 s frame-time. We applied the data process- 
ing techniques recommended by the Chandra X-ray 
Center (CXC) which include removal of hot pixels 
and non-X-ray events. The standard CXC software 
pipelines randomize the event positions by ± 0"246 
(1 ACIS pixel = 0"492), in each spatial coordinate. 
This randomization is done to avoid aliasing effects 
noticeable in observations with exposure times less 
than ~ 2 ks. Since the exposure time for the Chan- 
dra observation of Q0957-I-561 is considerably longer 
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than 2 ks, we reprocessed the data with the CXC 
software tool ACISPROCESSEVENTS without random- 
izing the event positions. This reprocessing led to 
an improvement of the 50 percent encircled energy 
radius from 0"A0 to 0".35. The CCD read-out streaks 
of the bright images A and B of Q0957+561 are re- 
moved in the spatial and spectral analysis. These 
streaks are produced from events that are recorded 
on the CCD during the transfer of charge into the 
CCD frame-store region. The background level was 
fairly constant throughout the observation with a 
mean value of 2.2 X 10~^ events arcsec"^ in the 
0.3-7 keV band. 

3. SPATIAL ANALYSIS; NEW X-RAY COMPONENTS 
IDENTIFIED 

3.1. The Lens Cluster 

In Figure 1 we show the Chandra image of the GL 
system Q0957-I-561. To enhance the presence of soft 
extended X-ray emission, the image was filtered to 
include only photons with energies ranging between 
0.5 and 3.0 keV. This choice of energy range is sup- 
ported by the analysis of the extracted spectrum of 
the cluster (see section 4.1) which indicates that the 
cluster spectrum is very soft with most of the de- 
tected counts lying below 3 keV. The image has been 
binned with a bin size of 0".5 and smoothed with the 
software tool CSMDDTH developed by Ebcling et al. 
(2000) and provided by the Chandra X-ray Center 
(CXC). CSMOOTH smoothes a two-dimensional image 
with a circular Gaussian kernel of varying radius. Ex- 
tended X-ray emission with a radius of ~ 30", en- 
compassing the quasar images, is clearly visible in 
Figure 1. The extended emission appears to be cen- 
tered slightly northeast of image B. This is seen if one 
takes the midpoint of the two outer contour levels of 
the X-ray emission shown in Figure 1. The midpoint 
is slightly to the cast of image B. 

To constrain the morphology of the soft extended 
emission we performed a 2D fit to the X-ray data 
using the model described below. The spatial model 
has the following three components: extended cluster 
emission, unresolved emission from the quasar images 
and background emission. 

(a) To describe the cluster brightness profile, we 
use a P model 

fir) = A[l + (^)2]-3/?+^ (1) 

where. 



x'2(l-e)2+y' 



x' = (x - xq) cos{6) + [y - yo) sin{9) 
y' = (y - yo) cos{9) - {x - xo) sin{e) , 

€ is the ellipticity of the model defined here as e = 
1 — q; q is the axis ratio; 9 is defined as the angle be- 
tween the major axis and west, and is measured west 
to north; xq, yo arc the positions of the cluster center; 
and (3 is the ratio of kinetic energy per unit mass in 
galaxies to kinetic energy per unit mass in gas. 

(b) To model the lensed quasar images, we use sim- 
ulated point spread functions (PSF's). The centroids 
of the model PSF's are fixed to the observed image 
centroids. The relative normalization of the PSFs 
is the flux ratio of the lensed images as determined 
from annular regions of inner and outer radii of 0"5 
and 2", respectively, centered on the images. These 
annuli were chosen to avoid the slightly piled-up cores 
of the images. The model PSFs appropriate for this 
observation were created employing the simulation 
tool MARX v3 (Wise et al. 1997) with an input spec- 
trum derived from the best-fit Chandra spectrum of 
the lensed images. Specifically, we used an absorbed 
power law with a column density of Nj^ = 0.82 x 10^° 
cm~2 (Dickey & Lockman, 1990) and photon indices 
of F = 2.08 and 1.94 for images A and B respectively 
(see section 4.2). Here-after the term best-fit is used 
to describe a result obtained by fitting a model to our 
data using an optimization technique to find the lo- 
cal fit-statistic minimum. The Levenberg-Marquardt 
optimization method is used for the spectral analysis, 
the Powell method is used for the spatial analysis of 
the cluster and the downhill simplex method is used 
by the fitting tool LYNX (section 4.2). 

(c) Finally in our model, we include a uniform 
background of 0.01 events per pixel obtained from 
a background region at the same distance from the 
aim point as the cluster. 

In determining the cluster properties, we omitted a 
10" radius region centered on the midpoint between 
the quasar images to avoid biasing the fit due to resid- 
uals in modeling the cores of the PSFs. For the an- 
nulus from 8" to 40" centered on the midpoint of the 
quasar images, we binned the image in 1" pixels and 
smoothed this with a Gaussian (a = 3") prior to per- 
forming the spatial fitting. The fits were performed 
with the CXC software package SHERPA. To evaluate 
the sensitivity of the results to the choice of assump- 
tions made for sizes of spatial windows and bin sizes 
of the data, we performed the spatial analysis over a 
wide range of input model parameters. Specifically, 
we binned the X-ray data with bin sizes varying be- 
tween 1" and 2". The annulus region used to model 
the cluster was varied with an inner radius ranging 
between 10" and 14" and an outer radius ranging be- 
tween 30" and 40". The errors quoted for the spatial 
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model parameters represent the maximum uncertain- 
ties obtained from all spatial fits. Results for our 
spatial fits arc presented in Tabic 1. Wc find the cen- 
ter of the mass of the lens cluster to be located at 
Aa = 4:".3t\i east, AS = 3'f5tJ| north of the core 
of image B. The fits indicate that the smoothed mass 
distribution of the cluster is close to spherical with an 
ellipticity e = 0.19lo:o|. The best-fit values for (3 and 
the core radius of the cluster are /3 = 0.47lQ Qg and 
ro = 15'J4l'^3 5 (62.5 kpc), respectively. Our values for 
f3 and ro, within the uncertainties, fall on the trend 
line of /? versus ro obtained previously for a suite of 
galaxy clusters (Jones & Forman 1999). 

3.2. The X-ray Jet of Image A 

We notice a faint feature in the Chandra image ex- 
tending northeast of image A that is aligned with the 
radio jet and appears to have a jet-like morphology. 
The X-ray jet is significant at the 4 a level. At a red- 
shift of 1.41 this is the most distant X-ray jet detected 
to date. An overlay between this X-ray jet feature 
and a 3.6 cm radio map of the jet in image A (radio 
data from Harvanek et al. 1997) is shown in Figure 2. 
The X-ray and radio jet morphologies appear to be 
somewhat similar. A comparison between the radio 
and X-ray jet profiles along the jet ridge line is shown 
in Figure 3. The X-ray jet extends NE about 8" from 
the core of image A with bright knots at 2"3, 4" and 
6" (here-after also referred to as knots A, B and C, 
respectively) from the core. Each of these knots coin- 
cides within 0'.'5 with radio knots at similar locations. 
The intensity of the X-ray knots appears to decay as 
a function of distance from the core in contrast to 
the radio knots which become brighter at larger dis- 
tances from the core. One possible interpretation of 
the former is that the jet flow is decelerating and 
the brightness change is due to aging of the higher 
energy electron population of the jet. The increase 
in the radio brightness along the jet may be due to 
an increase of the magnetic field strength. A similar 
anti-correlation between radio and X-ray profiles was 
recently reported in 3C273 (Sambruna et al. 2001; 
Marshah et al. 2001). Lens models for Q0957-F561 
predict a small magnification gradient along the jet. 
We estimate that the magnifications for the knots B 
and C are 1.9 and 1.6, respectively. Therefore, the 
magnification gradient does not greatly change the 
X-ray and radio brightness profiles. In section 4.3 we 
briefly discuss possible mechanisms that may explain 
the origin of the X-ray jet emission and present the 
spectral energy distribution for knots B and C. 

4. SPECTRAL ANALYSIS 
4.1. A Cool Lens Cluster 



The spectrum for the cluster of galaxies was ex- 
tracted from a 40"-radius circle centered on the X-ray 
determined cluster center. Beyond this radius, clus- 
ter emission is not detected above the background. 
Three-arcsecond-radius circles centered on images A 
and B were excluded. Emission from the X-ray jet 
was also excluded by omitting a rectangle 2".5 by 7". 
However, even beyond the 3" radius, mirror scatter- 
ing of the bright lensed images produces a significant 
contamination of the cluster spectrum, particularly 
at hard energies. We estimate the spectrum of the 
contamination produced by images A and B to the 
cluster spectrum by performing simulations with the 
raytrace simulator tool MARX provided by the CXC. 
We simulate two point sources centered at the lo- 
cations of images A and B with input spectra and 
normalizations derived from our spectral analysis of 
these images. We find a total of about ~ 200 counts 
due to quasar emission 0.8% of the total counts 
from images A and B) within the cluster extraction 
region. The CCD background and quasar contamina- 
tion are subtracted from the cluster region resulting 
in the spectrum shown in Figure 4. A net total of 
about ~ 600 X-ray events originate from the clus- 
ter. The data were fit with the spectral analysis tool 
XSPEC (Arnaud 1996). Wc model the cluster spec- 
trum with a Raymond - Smith thermal plasma mod- 
ified by absorption due to our Galaxy. The Galactic 
column density is fixed at the value of N// = 0.82 
X 10^*^ cm~^ for all spectral fits performed in our 
analysis. We assumed a typical value for the metal 
abundances in clusters of 30% solar (see, e.g., Hen- 
riksen 1985; Hughes et al. 1988; and Arnaud et al. 
1987). Due to the relatively low signal-to-noise ratio 
of the cluster spectrum the metal abundances of the 
cluster cannot be constrained within useful limits. If 
we let the abundance be a free parameter we obtain 
a best-fit value for the abundance of A = 0.19lQ^g 
(90% confidence level) and a best fit temperature of 
2.1keV. We note that the choice of metal abundance 
(between and 100% solar) has little effect on the 
temperature determinations. The abundance ratios 
used in the Raymond - Smith thermal plasma emis- 
sion model are those of Anders Sz Grevesse (1989). 
Our best-fit model is shown in Figure 4. The uncer- 
tainty in the calibration of ACIS S3 below 0.5 keV 
contributes to the large residuals between 0.4 and 
0.5 keV. We find a temperature for the cluster of Tg = 
2.09l|j keV at the 90% confidence level. The spec- 
tral fine features between 0.7 and 0.9 keV (observed 
frame) correspond to a redshiftcd z= 0.36 complex of 
Fe L lines. Our spectral analysis thus confirms that 
the extended emission originates from the lens. The 
relatively low cluster temperature of ^ 2 keV and 
our values of /? and the core radius from section 3.1, 
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within their uncertainties, are consistent with the ob- 
served correlations between temperature, (3 and core 
radius obtained previously for a large sample of clus- 
ters of galaxies (Jones & Forman 1999). We estimate 
that ~ 3% of the counts originating from cluster emis- 
sion are excluded from our extracted spectrum due to 
the regions used in our analysis. Correcting for this 
effect we find a 2 - 10 keV cluster luminosity of 4.7 
X 10^^ erg s~^. Our estimated values for Lx and Tg 
are consistent with the empirical Lx - Tg correlation 
between cluster luminosity and temperature of, Tj^^v 
= 0.13(+0.08,-0.07) X 2,g-35i±o o68)^ ^^lere L40 is the 
measured 0.5-4.5 keV luminosity in units of 10^^ erg 
s"! (Ho = 50 km s"^ Mpc"^) (Jones k Forman 1999; 
Markevitch, 1998). Specifically, the Lx - Tg correla- 
tion yields a value of Tg ~ 2 keV for our value of L40 
= 2.3 X 10^. 

4.2. The Spectra of Images A and B 

As mentioned in section 2, we expect the spectra 
of images A and B to be slightly piled- up. Wc used 
two independent methods to account for pile-up. 

In the first method we extracted spectra from an- 
nuli centered on each image with inner and outer 
radii of and 3", respectively. Pile- up is signifi- 
cantly reduced beyond the O'fS radius, in the wings 
of the PSFs. One of the drawbacks of the annulus 
method is that only ~ 40% of the total quasar counts 
are considered in the spectral fit, thus, leading to 
larger uncertainties in the values of the estimated pa- 
rameters. We corrected the ancillary files provided 
by the CXC to account for the energy dependence 
of X-ray scattering in the Chandra mirrors within 
the selected annuli. The energy dependent correc- 
tion function was evaluated by performing raytrace 
simulations with the software tool MARX. We modeled 
the spectra with power-laws modified by absorption 
from our Galaxy. The fits are statistically acceptable 
= 113(138) and x^(z^) = 109(105), for images 
A and B, where u is the number of degrees of free- 
dom) with best-fit values for the photon indices of 
images A and B of 2.06l[j:[]^ and 2.06l[]:[5f, respec- 
tively (90% confidence errors). The flux ratio B/A 
is 0.74 lb 0.02, consistent with the observed flux ra- 
tio in the radio band of 0.76 ± 0.03 (VLBI A 13cm, 
Falco et al. 1991) and 0.72 ± 0.04 (VLA A 6cm core, 
Conner et al. 1992). Previous measurements of the 
X-ray flux ratios of 0.3 it 0.1 and 1.5 it 0.2 made with 
EINSTEIN and ROSAT, respectively, (Chartas et al. 
1998) differed significantly from this value, suggesting 
the presence of microlensing for those epochs. 

For our second approach we used the forward fitting 
tool LYNX developed at PSU (Chartas et al. 2000). 
Spectra were extracted from circles centered on each 
image with radii of 3". We found the photon indices 



for A and B to be 2.08 ± 0.03 and 1.94 ± 0.03 (90% 
confidence errors). The 0.5-10 keV X-ray fluxes of 

images A and B corrected for pile-up are 11.6 it 0.9 
X 10"^2 erg s"^ cm-2 and 8.8 ± 0.7 x 10"^^ erg 
s~^ cm~^. The errors in the estimates of the photon 
indices in both methods do not include systematic 
errors due to the uncertainties in the ACIS detector 
quantum efficiency and energy rcsponce. We note 
that the "annulus" method uses detector response 
matrices provided by the CXC whereas LYNX links to 
a Monte Carlo simulator of ACIS developed by the 
ACIS instrument team (Townsley et al., in prepara- 
tion). We do not detect any line features in their 
spectra. Combining the spectra of images A and B, 
we place a 95% confidence upper limit of ~ 60 eV 
(observed-frame) on the equivalent width of an in- 
trinsically narrow fiuorescent Fe Ka line at 6.4 keV 
(2.66 keV observed-frame). 

4.3. The Spectral Energy Distribution of the Jet 

In Figure 5 wc show the spectral energy distribu- 
tion(SED) of the two brightest knots B and C of the 
jet. The radio fiux density, F^, at 3.6 cm for jet A 
was based on the value reported by Harvanek et al. 
(1997). The jet knots have not been detected in the 
optical band; we therefore placed upper limits on the 
optical flux density of the knots based on deep HST 
observations of this field performed by Bernstein & 
Fischer (1997). We find the Sa upper limit on vFi, 
at 5500 A to be 5.6 x 10"^^ erg s"^ cm'^. The 
X-ray spectrum was extracted from a 2.5" by 7" rect- 
angle. The background was extracted from several 
similar rectangles located at the same distance from 
image A but at different azimuths. A net total of 
50 ± 8 X-ray events were those ascribed to the jet. 
We modeled the spectrum with a power-law modified 
by absorption due to our Galaxy. We find a best-fit 
value for the photon index F of 1.87 ± 0.6 (90% con- 
fidence errors) including systematic errors estimated 
from using different background regions. We find I'Fu 
at 1 keV for knots B and C of 8.9 ± 2.8 x lO'^^ and 
2.7 ± 0.9 X 10-16 erg s'^ cm-^, respectively (90% 
confidence errors). A simple synchrotron model with 
a single power law distribution is not consistent with 
the SED of the jet in Q0957A. Given the high red- 
shift of the jet, the energy density of the Cosmic Mi- 
crowave Background (CMB) is enhanced by a factor 
of (1 -|- z)'^ = 33.7 over the local value. We there- 
fore modeled the spectra of knots B and C with a 
CMB model, in which cosmic microwave photons are 
inverse Compton scattered by synchrotron-emitting 
electrons in the jet. The solid and dashed lines in 
Figure 5 are the fits of the CMB model to the SEDs 
of knots B and C, respectively. We are able to fit 
both knots well with very moderate Doppler beaming 
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(see Figure 5). We conclude that a plausible mecha- 
nism to explain the X-ray emission is inverse Comp- 
ton scattering of CMB photons. Such a mechanism 
was shown to be consistent with the peculiar SED of 
the jet in PKS 0637-752 (Tavecchio et al, 2000). 

5. MASS DISTRIBUTION AND CONVERGENCE OF 
THE LENS CLUSTER 

For isothermal and spherical clusters of galaxies the 
equation for hydrostatic equilibrium can be solved for 
the virial mass, Mgrav{< r), within a radius r, 



Mgravi< r) = - 



kTr d In p 



(2) 

jxmpG d In r 

where k is Boltzmann's constant, /xm^ is the mean 
molecular weight of the cluster gas, and p{r) is the 
cluster gas density at r. 

Assuming a (3 model for the density profile of the 
hot gas (e.g., Jones & Forman 1984) we obtain the 
equation for the total mass of the lens within a radius 



Mgrav{< r) 



3pkT 



(3) 



p.mpG [1 + (^)2] ' 

By incorporating the best-fit spatial and spectral 

parameters from our present analysis (see Tables 1 
and 2) we find the total cluster mass within a radius 

of 1 /l75-^ MpC to be Mgrav = ^-^^H X 10^^ Mq. 

In Figure 6 we show the total cluster mass within a 
radius r as a function of radius. The shaded region 
indicates the allowed values for the cluster mass in- 
cluding the uncertainities obtained from the spatial 
and spectral fits to the lens cluster (see Tables 1 and 
2). 

These mass estimates were used to evaluate the 
convergence parameter k{x), 



KiX) 



S(x) 



where T,{x) is the surface mass density of the lens 
cluster as a function of the cylindrical radius x (Char- 
tas et al. 1998) and Scr is the critical surface mass 
density (see, e.g., Schneider, Ehlers & Falco 1992), In 
Figure 7 we plot the convergence parameter k{x) as 
a function of distance from the cluster center. The 
thick solid line corresponds to the best-fit spatial and 
spectral parameters. The largest contributor to the 
uncertainty in our estimate of k{x) is the weak con- 
straint on the temperature of the cluster. To illus- 
trate this we have plotted the uncertainty in k(x) as- 
suming 68% (dotted lines) and 90% (dashed lines) 
confidence intervals for the temperature. We also 
chose cluster limits ranging from O.Trsoo and 1.4r5oo, 
where r5oo is the radius in which the mean over- 
density is 500, and rsoo = 1.58(Tx/10 keV)^/^ f^-i 



Mpc ~ 0.71 h^^ Mpc (Mohr, Mathiesen, k Evrard, 
1999). The uncertainity in k introduced by assuming 
the cluster limit to range between 0.7 and 1.4r5oo is 
only ~ 2%. Our assumption of spherical symmetry of 
the mass distribution has little effect on the estimate 
of the convergence parameter. Calculations of the 
mass distribution of the nonsperical cluster of galax- 
ies A2256 (axis ratio ~ 1.6) by Fabricant et al. (1984) 
showed that radially integrated mass estimates were 
negligibly affected by including an oblate or prolate 
cluster geometry. At the best-fit location of images A 
and B with respect to the cluster center, we estimate 
the convergences assuming the 90% confidence range 
in cluster temperature and find ka = 0.221qq7 and 
Kb = 0.211q respectively. Our estimated range of 
uncertainty in k does not include possible systematic 
errors arising from the estimation of the the surface 
mass density of the galaxy cluster through the use of a 
hydrostatic, isothermal (3 model. Evrard, Metzler & 
Navarro (1999) have performed simulations to inves- 
tigate the accuracy of galaxy cluster mass estimates 
based on X-ray observations. They found that esti- 
mates based on the hydrostatic, isothermal (3 model 
are unbiased estimates of the mass of relaxed clusters 
with standard deviations of less than 30%. Recent ob- 
servations of clusters of galaxies with XMM-Newton 
indicate that the cluster temperature profiles are re- 
markably isothermal beyond the central cooling flow 
(e.g., Arnaud et al. 2001). We therefore expect that 
our estimated uncertainity in k is dominated by the 
large uncertainity in the estimated value of the clus- 
ter temperature and not systematic errors introduced 
from the use of the isothermal (3 model. 



6. DISCUSSION AND CONCLUDIONS 

Detailed lens models of Q0957-I-561 include the 
cluster's contribution to the lensing potential by ex- 
panding the potential of the cluster in a Taylor series 
and keeping terms of up to third order (Kochanck 
1991, Bernstein k Fischer 1999, Kecton et al. 2000). 
The second order term of the expansion represents 
the shear from the cluster and can be expressed as 
7 = (i+/3^'^^)3/2 where, Prd = rc/dc, Vc is the cluster 
core radius and dc is the distance from the cluster to 
the galaxy Gl (see, Kochanck, 1991). Recent models 
of Q0957-I-561 that incorporate additional constraints 
from observed arcs produced by the lensing of the 
quasar host galaxy imply a cluster shear amplitude, 
7, that is relatively small compared to the conver- 
gence parameter k (Keeton et al. 2000). Other ob- 
servations of Q0957-I-561 also favor values of k >> 7 
(Fischer et al. 1997; Romanowsky & Kochanek 1999; 
Chartas et al. 1998). As pointed out by Keeton et al. 
(2000), one way of producing 7 << k is to have the 
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distance from the center of the cluster to the center 
of the lens galaxy be: less than the core radius for a 
cluster with a singular isothermal ellipsoid potential, 
or: less than the "universal" scale length for a clus- 
ter potential that follows the dark matter profile of 
Navarro, Frenk, and White (1996). Measurements at 
that time, ie., the late 1990s, indicated that the clus- 
ter center was located several core radii away from 
the lens (Fischer et al. 1997; Angonin-Willaimc et 
al. 1994) and therefore were not consistent with the 
large value of n relative to the amplitude of the shear. 

The Chandra observation of the cluster resolves 
this apparent discrepancy. We find the distance be- 
tween the center of mass of the cluster and the cen- 
ter of the lens galaxy Gl to be dc = 4'f8 consider- 
ably smaller than suggested from previous estimates. 
The best-fit position angle (North through East) and 
core radius of the cluster are ~ 59° and 15'f4l3|, re- 
spectively. The Chandra observations of Q0957+561 
therefore indicate that the cluster center is located 
within a core radius of the lens galaxy Gl. Using the 
values for Tc, d^, and k provided by our analysis of the 
Chandra observation of Q0957-I-561 we find that the 
cluster shear amplitude is 7 = 0.01 ± 0.009 (90% con- 
fidence level) , consistent with the recent model results 
of Keeton et al. (2000). The Ciiandra observation of 
Q0957+561 has eliminated several uncertainties in- 
troduced in our previous analysis of the ROSAT HRI 
data of this system. Specifically, the temperature, the 
core radius, and cluster shape are now determined 
more reliably. Previous estimates of k are unreliable 
due to the large distances between the center of the 
cluster and the galaxy Gl assumed in these analyses. 
To evaluate the implication of our estimate of k on 
the Hubble constant, we write Hq = Hq{1 -k) = 100/t 
km s~^ Mpc~^ , where we define the average conver- 
gence parameter of the cluster at the image locations 
as K = {ka + k_b)/2. Due to the mass-sheet degen- 
eracy problem, lens models provide only the model- 
dependent value h/{l - k). Our present observation 
of the lens cluster results in a ~ it 14% uncertainty in 
the quantity (1 - k). The resulting uncertainty in Hq 
is ~ ± 29% , where we have assumed an uncertainty 
of ~ lb 25% in Hq based on the analysis of Keeton et 
al. (2000). We anticipate that the next generation of 
lens models for Q0957-I-561 may provide tighter con- 
straints on Hq when the constraints on the location, 
ellipticity and convergence of the cluster, based on 
the Chandra observations, are incorporated. 

As described in section 5 the largest contributor to 
the present uncertainty in k obtained by the X-ray 
method is the temperature of the intracluster gas. 
Given that the cluster is located near the quasar im- 
ages, we may optimize a future observation of Q0957 
by moving the telescope aim point closer to the clus- 



ter center to improve the spatial resolution and re- 
duce the contamination of the cluster by the bright 
images. This change will reduce the uncertainty of 
the estimates of the spectral and spatial parameter 
values for the cluster. 

We would like to thank Daniel Harris and William 
For man for helpful discussions. We acknowledge fi- 
nancial support by NASA grant NAS 8-38252 and 
support from the Smithsonian Institution. 
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TABLE 1 

Results Prom Spatial Fits To Lens Cluster 





yo 






n 


n 




n 




o r.+l.2 


0.471°:°^ 


15.4^1:1 0.19i°:°^ 



NOTES- 

The probability distributions for the best-fit model 
parameters yo and e derived from our error analysis 
are not Gaussian and any value of yo and e within the 
quoted errors should be considered to have similar 
likelihood. The errors quoted for the spatial model 
parameters, represent the maximum range of the pa- 
rameter values obtained from the suite of spatial fits 
performed in our sensitivity analysis. 
" xq and yo arc the separations Act (a positive value 
indicates east of B) and AS (a positive value indicates 
north of B) of the center of the lens cluster from the 
core of image B. 

^ /3 is the ratio of kinetic energy per unit mass in 
galaxies to kinetic energy per unit mass in the gas. 
P is determined from the fit to the surface brightness 
of the lens cluster. 

ro is the core radius of the lens cluster. 

e is the ellipticity of the cluster defined as e = 1 - 
q, where q is the axis ratio. 



9 



TABLE 2 

Values of Model Parameters Determined from fits to the Spectrum of the Lens Cluster 



Model 'Nh 
cm 


Te 

keV 


Fx 

erg s~-^ cm~^ 


erg s~-^ 




ABS + RS 0.82 X lO^o 


2.09l°jt 


l.lx 10-1^ 


4.7 X 10^2 


1.3(43) 



NOTES- 

The spectra are described by a thermal Raymond - 
Smith model plus absorption due to cold material at 
solar abundances fixed to the Galactic value. The 
absorbed flux Fx and unabsorbed luminosity Lx of 
the lens cluster are estimated for the energy range 
between 2 and 10 keV. The reduced chi-squared is 
defined as xi = X^/^j where the number of de- 
grees of freedom, is given in parentheses. 
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Fig. 1. — An adaptively smoothed X-ray image of the gravitational lens system Q0957-I-561. To improve the signal- 
to-noise of the lens cluster we selected X-rays with energies ranging between 0.5 and 3 keV. The lens cluster is resolved 
and appears to be centered slightly northeast of image B. The contour levels are 0.1, 3.5x 10"'*, 3.5x 10~^, 2.2x 10~^, 
and 1.3x 10~^ of the peak emission of image A. The white circles indicate the optical locations of the galaxies in the 
cluster at z = 0.36. Galaxy positions are taken from Angonin-Willaime et al. (1994). North is up east is left. 
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VLA and Chandra - Q0957 JET 




Right Ascension (J2000) 



Fig. 2. — An overlay of the VLA and Chandra images of Q0957+561. The contours represent the radio 3.6 cm VLA 
image of Q0957+561 provided courtesy of Harvanek et al. (1997). The X-ray and radio data are binned with a bin size 
of 0"048 on a side. For presentation purposes the X-ray data are smoothed with a circular Gaussian of radius a = 0".25 
resulting in an effective resolution of cr 0"33. The overlay clearly shows that the shape and angular structure of the 
X-ray and radio jets are similar. The radio contour levels are 0.11, 3.5x 10~^, 1.1 x 10~^ and 3.5x 10~^ of the peak 
emission of image A. The beam width is 0'.'24 x 0"17. The X-ray image contains photons with energies ranging between 
0.5 and 10 keV. The brightness profiles of the X-ray and radio jets differ significantly as shown more clearly in Figure 3. 
North is up east is left. 
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Fig. 3. — The radio and X-ray intensity profiles are plotted as a function of angular distance from the core of image A 
along the jet ridge line. The X-ray profile provides counts in 0'.'288 increments integrated ± 0"48 perpendicular to the 
jet ridge line and the radio profile provides the 3.6 cm fiux density in 0"096 increments integrated ± 0"48 perpendicular 
to the jet ridge line. Radio counterparts within 0'.'5 of the X-ray bright knots located 4" and 6" from the core are clearly 
visible. Whereas the angidar structure of the X-ray and radio jets are similar, the intensity distributions along the jet 
appear to be anti-correlated. In particular, the X-ray brightness decreases along the jet, whereas the radio flux density 
of the knots increases with distance from the core (see text for interpretation). 
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Fig. 4. — (Top panel) The X-ray spectrum of the lens cluster along with the best-fit thermal Raymond-Smith model 
(see Table 2). The model is consistent with a plasma temperature of ~ 2.1 keV. The detected spectral feature between 
0.7 and 0.9 keV (observed-frame) correspond to a redshifted z = 0.36 complex of Fe L lines, confirming the detection of 
the lens cluster. The uncertainty in the calibration of ACIS S3 below 0.5 keV contributes to the large residuals between 
0.4 and 0.5 keV. (Lower panel) Residuals in units of standard deviations with error bars of size la. 
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Fig. 5. — The spectral energy distributions (SEDs) of jet knots B (open circles) and C (filled circles). The radio 
flux densities were extracted from the radio A3. 6 cm VLA data kindly provided by Michael Harvanek of Apache Point 
Observatory. The optical upper limit (cross) is based on published results from a deep HST observation of the field 
(Bernstein & Fischer 1997). The values of uFi, of the jet knots in the X-ray band are at least an order of magnitude 
above an extrapolation of the radio and optical data points, siiggesting that the emission in the radio, optical and X-ray 
bands is not consistent with a single synchrotron model with a single power law distribution. The solid and dashed lines 
represent fits of the CMB model (see text) to the spectra of knots B and C, respectively. The CMB model, in which 
cosmic microwave photons are inverse Compton scattered by synchrotron-emitting electrons in the jet, has the following 
parameters. Knot B: The emission region is assumed to be spherical with a radius of i? = 1 x 10^^ cm, a magnetic field 
intensity of _B = 6 x 10^'"' G, and a Doppler factor of 5 = 1.4. The electron distribution is assumed to be a power-law 
with extremes of jmin = 10 and jmax = 3 x 10^, a slope of n = 2.6, and a normalization of K = 6.4 x 10~^ cm~^. Knot 
C: The CMB model for knot C has the same parameters as knot B except K = 4 x 10~^ cm~^ and d = 1.05. 
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Fig. 6. — The total mass of the lens cluster within a radius r. The dashed line corresponds to the best-fit spatial and 
spectral parameters. The shaded region indicates the allowable range of the cluster mass for the estimated uncertainties 
of the best-fit spatial and spectral parameters. 
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Fig. 7. The convergence parameter k of the lens cluster as a function of distance from the chister center. The thick 
sohd hne corresponds to the best-fit spatial and spectral parameters. The largest contributor to the uncertainty in the 
present measurement of k{x) is the weak constraint on the temperature of the cluster. To illustrate this weakness we 
have plotted the uncertainty in k{x) assuming 68% (dotted lines) and 90% (dashed lines) confidence intervals for the 
temperature. We also chose cluster limits ranging from 0.7r5oo and 1.4r5oo, where r^oo is the radius in which the mean 
over-density is 500, and rsoo = 1-58 hjr^ Mpc (Tx/lO keV)!/^ _ g.Tl h^^ Mpc (Mohr, Mathiesen, & Evrard, 1999). The 
solid vertical lines indicate the distances of images A and B from the cluster center. 



